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We propose the use of a superconducting charge qubit capacitively coupled to two resonant nanomechanical 
resonators to generate Yurke-Stoler states, i.e. quantum superpositions of pairs of distinguishable coherent states 
180° out of phase with each other. This is achieved by effectively implementing Kerr nonlinearities induced 
through application of a strong external driving field in one of the resonators. A simple study of the effect 
of dissipation on our scheme is also presented, and lower bounds of fidelity and purity of the generated state 
are calculated. Our procedure to implement a Kerr nonlinearity in this system may be used for high precision 
measurements in nanomechanical resonators. 
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I. INTRODUCTION 

Quantum nonlinear dynamics is an important topic in 
physics. In quantum optics, nonlinear interactions have been 
widely used to generate nonclassical field states, such as 
squeezed or sub-Poissonian light 1 . A special class of opti- 
cal nonlinearity results in an intensity dependent phase shift, 
commonly known as the Kerr effect. In the single mode case, 
the time evolution of an initial coherent state, under the influ- 
ence of such a Kerr medium and very low loss, will evolve 
into a quantum superposition of two coherent states 180° out 
of phase with each other. This was first discovered by Yurke 
and Stoler 2 , and since then such states have been called Yurke- 
Stoler states. A single-mode Kerr medium preserves the pho- 
ton statistics but modifies the quadrature uncertainties gener- 
ally leading to squeezing^. 

There is great interest in observing this quantum nonlin- 
ear couplings in solid state systems. This would allows us 
to deepen our current understanding of the classical-quantum 
frontier by studying how long can superpositions of meso- 
scopically distinct states survive in such systems. Some in- 
teresting proposals involving nanomechanical resonators have 
been published during the last years. In one scheme 3 , the use 
of a time dependent drive in a Cooper pair box coupled to a 
nanomechanical resonator is shown to generate a number of 
nonlinear Hamiltonians for the latter. By parametrically driv- 
ing a nanomechanical resonator capacitively coupled to a su- 
perconducting coplanar waveguide one can generate interest- 
ing nonlinear Hamiltonians suitable for generation and detec- 
tion of squeezed states as proposed in 4 .. In this system, entan- 
gled states in temperatures up to tens of milliKelvin may be 
achieved as discussed irA Nanomechanical oscillators have 
also been shown to be feasible for coupling to other impor- 
tant physical systems besides Cooper pair boxes or microwave 
fields of coplanar wave guides. Nanomechanical resonators 
may, for instance, be coupled to Bose-Einstein condensates 6 -, 



trapped ionsiS. or spin degrees of freedom of a sample of neu- 
tral atoms in the gas phased 



In this paper, we propose a theoretical scheme to engineer 
Kerr Hamiltonians using a system composed of a Cooper pah- 
box capacitively coupled to two resonant nanomechanical res- 
onators. We show in Sec. |II]that such nonlinear Hamiltonians 
can be achieved in a dispersive regime by appropriately choos- 
ing the system's parameters and by using a properly tuned 
strong classical field in one of the resonators. The integration 
of superconducting qubits with nanoresonators is an impor- 
tant topic and has been previously considered i n 10 ' 11,12 ' 13 ! 14,15 . 
We start from a well known Hamiltonian describing the in- 
teraction between a charge qubit and two resonant nanome- 
chanical resonators in a quantum regim e) 14 ' 16 ' 17,18 , and we then 
include an external driving in one of the oscillators. By con- 
sidering the regime of intense driving, we show that a non- 
linear Kerr-type effective Hamiltonian may be obtained. This 
Hamiltonian is induced by the common coupling of the res- 
onators with the qubit and intense external driving. This is the 
central result of this paper, and as an application, we show in 
Sec. [HI] how to generate the Yurke-Stoler state in the normal 
modes of the nanomechanical resonators. We also discuss the 
zero temperature decoherence in a particular regime of relax- 
ation, and evaluate both the fidelity and the purity of the gen- 
erated superposition state. Finally, we would like to point out 
that the ability to implement Kerr nonlinearities in nanome- 
chanical resonators has recently been shown to find applica- 
tions also in high precision measurements. In a recent paper 19 , 
Woolley et al. have proposed a new protocol for high preci- 
sion measurement in a nanomechanical resonator that makes 
explicit use of such nonlinearities. This might be a potential 
application for the results presented in this paper. In Sec. [IV] 
we draw some conclusions. 
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II. THE MODEL AND THE KERR TYPE INTERACTIONS of the qubit to new variables a a -> a c , 



The simplest charge qubit, the Cooper pair box (CPB), con- 
sists of a small superconducting island with an excess number, 
n, of Cooper-pairs, connected by a tunnel junction (capaci- 
tance Cj and Josephson coupling Ej) to a superconducting 
electrode. External control is achieved by the application of a 
voltage gate V g coupled to the CPB via a gate capacitor with 
capacitance C g . More details can be found in the review^. 
For specific qubit proposals and decoherence analysis see^l. 
In our study, we will assume that the CPB is coupled ca- 
pacitively to two nano-eletromechanical systems (NEMS) 18 , 
as depicted in Fig.(Q~|i. In the two level approximation for 
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FIG. 1 : Two nanomechanical resonators (with lowering operators a 
and b) are capacitively coupled to a Cooper pair box (CPB). One of 
the oscillators is driven by a classical force. 

the CPB, the capacitive coupling between the qubit and two 
NEMS is described by the Hamiltonian^ 



H = fkua^ a + hujb^b -\ — ^-a z -\ — — a x + hXi(a + a})a z 



+HX 2 (b + tf)a z + hg{ae w ^ + oV*"'*), 



(1) 



where a, aS are the raising and lowering operators for the 
driven NEMS^ 2 -, b, are the raising and lowering operators 
for the other NEMS with the same resonance frequency ui, 
and g represents the amplitude of the external nanomechani- 
cal drive (frequency u> e ). The parameters appearing in ([TJ are 
given by 



Csdi V 2muj ' 



huj = -4E c (l-2n ff ), (2) 
HA = -2£. 7 cos(7r0/0 o ), (3) 

(4) 

where G g .i is the capacitance between the CPB and i-th 
nanomechanical bias gate, Cs is the total capacitance, di is 
the distance between the i-th nanomechanical bias gate and 
the CPB, and m is the mass of the NEMS. The couplings can 
be made different by varying, for example, the distances di or 
applying DC voltages to the resonators 2 ^. Our goal now is to 
show how the application of the external driving field may be 
used to engineer nonclassical states. We first make a rotation 



Now, moving to a rotating frame with frequency u> e , and set- 
ting Cl = u! c and 5 = u> — u e , we get 

H = hSata + hStfb + hgia + a!) + hX 1 (ae- iw ° t + 
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HX 2 (be- l ^ t + b^e luJ " t )[cos9a z - sm9(a+e lul ' t 
+<r_e- iw '*)]. (9) 

We can now make the rotating wave approximation, to get the 
interaction picture Hamiltonian, 
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An interesting situation appears when one takes the dispersive 
approximation (| S | 3> A i, A2, g) for the above Hamiltonian 
(applying similar methods to those described in2i). In this 
regime, the Hamiltonian ( fTOb may be approximated by 



H = hQa^aa z 
where 
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The Hamiltonian (fTTT i can be diagonalized by using new 
bosonic composite operators a\ = (cos + sin ^b) and 
a 2 = (— sin T^a + cos ^b) with appropriate choice for 7. We 
set from now on il = x, i.e. Ai = ±A2, since for this case 
the simple choice 7 = § solves the problem. In terms of the 
new operators a\ = 2~ 1 / 2 (a + b) and a 2 = 2~ 1,/2 (a — 6), H 
is written (setting h = 1) as 
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where £ = — j 1 sin 8. We will now show that in the regime 
I A I ^S> |C|> a Kerr type Hamiltonian can be generated. From 
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2 2 

( = — sin 6, we see that | A| ^> |£| implies that we must 
have g ^> 2AiSin(9 (i = lor 2), i.e a strong driving (g 
Xi). To make this clear, lets us assume Ai = A2 and A > 0. 
By transforming H + to an interaction picture with respect to 
Aov, one obtains 



V + (t) = ^{a\ ai [(a z 



ia y )e + [a z + ia v )e J}. 



(17) 



Now, if one defines the operator A = a\ai (a z — ia y ) and the 
constant A = |, the above Hamiltonian will read V+(t) = 
\(Ae l2At + At e ^ 2At ). It can be showr^i that for A » A, 
the effective Hamiltonian V_£ ff = h^x [A 1 A^] can be used. By 
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evaluating this commutator, one finds 
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where /.t = ( 2 /2A. Remarkably, this Hamiltonian mimics the 
single mode Kerr effect. If the CPB is prepared in an eigen- 
state of a x , the bosonic mode will follow a decoupled evolu- 
tion under the nonlinear Hamiltonian /j,(a[ai) 2 . Going back 
to the definitions, one can see that the magnitude of the non- 
linearity [i is in fact controlled by the system parameters Ai 
(coupling constant for the interaction of resonator a with the 
qubit), g (related to the amplitude of the classical driving), and 
S (detuning between driving field and nanoresonators). Thus, 
it is possible to control the the intensity of the present Kerr 
type effect, which is always important in the applications. 



III. YURKE-STOLER STATE AND INCLUSION OF 
DISSIPATION IN THE NEMS 

Consider now the initial preparation, = 
|a)a|o ; )fc|+) a ;, i.e both resonators in coherent states with 
the same amplitude a, and the CBP in an eigenstate of a x 
with eigenvalue equal to one. In the transformed space of 
the composite modes a\ and a2, this initial state becomes 
|V>(0)) = |v / 2a)i|0)2|+)a;- It means that the composite 
mode-1 is initially in a coherent state \ct\ = s/2a)i, mode-2 
in the vacuum state \ct2 = 0)2, and the qubit in the eigenstate 
of a x corresponding to the eigenvalue 1. For this initial 
condition, Hamiltonian ( fT8l leads to the following time 
evolved state: 
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with ol\ = \/2a. For an interaction time tj such that /itj = 
tt/2, the state (Q2]l evolves to \4>i(ti)) = |YS)i|0) 2 |+), 
where 



|YS>, 



|«i)i +i\ - 
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(20) 



is the Yurke-Stoler state. We remark that no measure- 
ment whatsoever was needed to generate this state, so 



this scheme is deterministic. If initially one prepares 
IV'(O)) = \ a )a\ — a )b\+)x an d choose Ai = — A 2 , a Yurke- 
Stoler state is generated in the mode-2. Many applications for 
superpositions of coherent states have been suggested in the 
quantum optics and quantum information literature 25 , along 
with a considerable variety of generation protocols^. 



Since we have performed a perturbation approach of the 
problem (effective Hamiltonians), it is now important to make 
a brief discussion about the experimental values of the param- 
eters and the feasibility of the regimes we used. From (O to 
( TTOt . we have realized a rotating wave approximation, and this 
is justified when Ai, A2 <C cj, f2, g. According to experimen- 
tal reference^, it is currently possible to achieve ui/2ir = 1.0 
GHz. For charge qubits, ordinary values for f2 are also about 
a few gigahertz 2 ^. In principle, the external driving g may 
also be of the same order or even stronger than ui and A. We 
have also demanded f2 = ui e , and this means that the fre- 
quency of the drive field is also of a few gigahertz. Taking all 
these into account, we see that the coupling constants Ai and 
A2 must be at most around a few megahertz for the rotating 
wave approximation to be valid. This seems not be a prob- 
lem since such coupling constants may be tuned by chang- 
ing the distance between the CPB and the nanoresonators or 
through additional DC voltages on the resonator. When going 
from ( [Tol l to ( fTTI ). we took the dispersive regime that demands 
Ai, A2 <IC S, Again, this might not be a problem since Aj de- 
pends on di. Finally, our last approximation corresponds to 
the regime of strong driving Ai, A2 -C g- This seems to be 
easy to achieve since g is externally controlled via a driving 
gate and do not depend on the fabrication features of the CPB 
or the resonators. 



It is well known that superposition states of this kind are 
easily corrupted in noisy or dissipative environment. For this 
reason, it is important to find a way to evaluate, at least ap- 
proximately, how our generation protocol is affected by such 
irreversible effects. A complete treatment of the problem 
would involve modeling the qubit decoherence and relaxation 
as well as different dissipative effects in the nanomechanical 
resonators. It is not our intention here to account for all these 
noise mechanisms. Instead, we will present one simple situ- 
ation which allows of a very illustrative exact solution. We 
consider the case in which both NEMS (with lowering opera- 
tors a and b) lose energy to their surrounding with decay rates 
K a and Kb, respectively. For simplicity, we will not include 
the qubit decoherence and relaxation. This is justified if the 
qubit decoherence times are longer compared to the resonators 
ones. At present, the charge qubits are notably more robust 
against decoherence and relaxation than the nanomechanical 
resonators. In this situation, and considering Ai = A2, the 
system master equation at zero temperature, when expressed 
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in terms of the mode operators, is written as 

®f>I -r / t \2 ~ I 

— = -i[fj,{a[ai) a x ,pi\ 

H 1 — (2aip/a{ - ajaipj - pia[ai) 

K a + Kb (r. ~ t t " - - t ^ 

H ^ [2a 2 pia 2 — a 2 a 2 pi — pia' 2 a 2 ) — 

^(i r/t t\ / t t\i 

^ K a i a 2 + a 2 ai)pi + pi(a[a 2 + a 2 a\)\ 

+^^[2(^ + ^)^(4 + 4)]. (2D 

We can see that the master equation contains extra terms 
due to the transformation to normal modes. The full treat- 
ment for arbitrary K a and Kb makes analytical progress quite 
difficult 2 ^ and a numerical calculation may be presented else- 
where. However, the simple regime in which both resonators 
decay with similar rates can be readily investigated. In this 
case, |/c + Kb\ \K a — Kb\, and we can drop the terms 
proportional to (/c — Kb). This assumption is realistic here 
since both resonators are assumed to be identical (same mass 
and natural frequencies). Even when this is not exactly the 



case, the quantities calculated below under the assumption of 
(k„ w Kb) will, at least, serve as an upper bound to the case 
in which the dissipation rates are disparate. 

We calculate here the degree of purity and the fidelity of 
the state, generated in such a noisy environment, as compared 
to the Yurke-Stoler state obtained in the ideal unitary case. 
Therefore, we need to take the same initial preparation used 
in the ideal case i.e. |^(0)) = |a) a |«)b|+).T- As this implies 
that mode-2 will be in the vacuum state, we need to consider 
only the terms in (fSTJ that contain operators for mode-1. The 
master equation ( f2TT > reduces to, 

^EL = -i[^(4ai) 2 ,/5/] + K,(2aip~ia\ - a\dxpi 

-pia\ai), (22) 

where k — (k q + k&)/4. An exact solution for ( 1221 using the 
Q-function approach is presented in 2 ^, but we will use a recent 
solution obtained directly for the density operator— which al- 
lows us to readily obtain the purity P = Tr[p 2 ], and fidelity 
F = (YS|p|YS). According to22, the solution of ([22} is 



Pl(t) = { ^1 Pn+k,m+k(0)e' 
fc,n,m— 

®|0> 22 (0| ® \+)xx(+\, 
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(23) 



where p n ,m(0) are the (Fock) matrix elements of the initial 
density matrix of mode-1. In Fig.©, we show the decay of 
fidelity (solid) and purity (dotted), as a function of the dimen- 
sionless parameter T = k/ p for a = 2, at the time for which 
the Yurke-Stoler state arises, i.e., pti = it/2. We can see 
that the fidelity is quite high (F > 0.99) for T < 10~ 3 . As 
expected, Figf2] reveals that the purity is more affected with 
increasing T than is the fidelity. However, it also presents sat- 
isfactory values for V w 10 -3 (P « 0.99). For more realistic 
values such as T « 10" 2 , we find F > 0.95 and P > 0.90. 

Finally, a few words about detection of superpositions of 
pairs of distinguishable coherent states is in order. This is an 
important topic, and several methods for detecting these states 
have already been proposed in the literatur e) 31 - 32 ' 33 . Among 
them, it seems that the most suitable method for the system 
treated here is the one presented in^ 2 -, whereby motional states 
of a single trapped ion have been experimentally determined. 
This method relies upon implementation of displacement op- 
erators and Jaynes-Cummings interactions to determine both 
the density matrix in the number state basis and the Wigner 
function. Thi ion techniques could be an alternative to detect 
the Yurke-Stoler state proposed in this paper, but it should be 
remarked that a CPB coupled to two NEMS has not yet been 
operated in strong coupling regime. 



IV. CONCLUSIONS 

To summarize, we have proposed a theoretical scheme to 
engineer a nonlinear Kerr Hamiltonian using superconducting 
charge qubits and nanoresonators. We have shown how such 
systems may be used to mimic a Kerr Hamiltonian. The for- 
mation of the Yurke-Stoler states in the composite mode of 
both resonators occurs naturally at an appropriate interaction 
time without needing to make a measurement on the system. 
For the case in which both resonators have equal decay rates, 
a simple exact expression for the total density matrix was de- 
rived. The present treatment, while not complete (more com- 
plex models of dissipation and noise could be considered), 
serves as an upper bound for the case in which the qubit deco- 
herence can be neglected. In this context, we have shown that 
the fidelity of the generated state can high for moderate values 
of the decay constants. As a final remark, recently Woolley et 
al^ have proposed a new protocol for high precision mea- 
surement in a nanomechanical resonator that makes explicit 
use of a Kerr nonlinearity. The method of the present pa- 
per could enable the use of linear nanomechanical resonators 
for such measurements instead of the intrinsically nonlinear 
nanomechanical resonators assumed iniS. 
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FIG. 2: Fidelity (solid) and purity (dotted) in function of the dimen- 
sionless parameter T = k/(j, for a = 2, at the moment the Yurke- 
Stoler state would be perfectly generated in the ideal lossless case. 
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